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Abstract. We give a number of algorithms for constructing unitary 
matrices and tight frames with specialized properties. These were pro- 
duced at the request of researchers at the Frame Research Center 
(www.framerc.org) to help with their research on fusion frames, the 
Kadison-Singer Problem and equiangular tight frames. 

1. Introduction 

We will give a sequence of algorithms for constructing unitary matrices 
(i.e. orthonormal bases) and tight frames with specialized properties. For 
unitary matrices, we just need to construct square matrices with rows and 
columns orthogonal and row and column numbers square summing to one. 
For tight frames, we need to produce M x N matrices with M > N which 
have orthogonal columns and whose column numbers square sum to a fixed 
constant. Then the row vectors of this matrix will form a tight frame for 
Hat. We will be interested in a number of special properties such as sparse 
orthonormal bases for hyperplanes, matrices with first row a constant, first 
two rows a constant, matrices made from just 2 constants, matrices built 
out of unitary matrices, matrices with a large portion of the norms of their 
vectors allocated to special positions in the matrix etc. We will work specif- 
ically with 2-tight frames since it is known that paving these is equivalent 
to solving the famous Kadison-Singer Problem [U|4j[5]. 

Recall that a family of vectors {fi}fLi ls a frame for a Hilbert space M.^ 
if there are constants < A < B < oo, called the lower and upper frame 
bounds respectively so that 

M 

A\\f\\ 2 <Y,\(f,fi)\ 2 <B\\f\\ 2 . 
i=l 

The frame is A-tight if A = B and Parseval if A = B = 1. It is equal norm 
of = c, for all i = 1,2,..., M, and unit norm if c = 1. The analysis 
operator of the frame is T : Hjv — > given by 

M 

i=i 
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where {e{}f£ l is the natural orthonormal basis of . The adjoint of the 
analysis operator is the synthesis operator which satisfies 



M 



T* ({a^ 1 )=E a ^- 



i=i 



The frame operator is the positive, self-adjoint invertible operator given by 



a; 



S=: T*Tf = J2(fJi)U 



i=l 



The synthesis operator of the frame is the matrix with the frame vectors as 
column vectors: 



A 



fi 



h 



M 



2. Sparse Orthonormal Bases for Hyperplanes 

For their work on sparsity the researchers at the Frame Research Center 
needed the sparsest orthonormal bases for hyperplanes in M. N which do not 
contain the span of any subset of the orthonormal basis. We believe that 
this class of hyperplanes is given by the following cases - although at this 
time we do not have a proof that these are the sparsest. 



dimension n=2: 



[1 1] 



dimension n=3: 



2 1 
-1 2 



dimension n=4: 



1111 

A 4 = 1 -1 
1-1 

The above has sparsity 8 and it is easily checked that any orthonormal basis 
for this hyperplane in M 4 has sparsity > 8. 

dimension n=5: Here we combine the earlier cases of A2 and A3 and then 
add a non-zero row on top. 



An 



The above has sparsity 11. 



* * * 
0**00 
** 
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dimension n=6: Here, we will combine the cases 2n = 2 with 2n = 4 and 
then insert a new top row. 



A 6 = 



1 


1 


H- 1 


1 


1 


1 


1 


1 


-1 


-1 








1 


-1 




















1 


-1 




















1 


-1 



The above has sparsity 16. 

dimension n = 7: Here we combine the earlier cases of n = 4 and n = 3. 



A 7 = 



* 


* 


* 


* 


* 


* 


* 


* 


* 


* 


* 











* 


* 























* 


* 























* 


* 


* 

















* 


* 



This has sparsity 20. 

dimension n=8: In this case we combine two of the cases 2n = 4 and add 
a new top row. 



^8 = 



h- 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 














1 


-1 


























1 


-1 


























1 


1 


-1 


-1 














1 


-1 


























1 


-1 



This has sparsity 24. 

The general cases here are: 

dimension 2n with n even: We combine two of the cases Am and add a 
new top row. 

. 2n ones 2n ones 

This has sparsity: 

2 [sparsity of ^2n] + 4n. 

dimension 2(n+l), n+1 odd: We combine the previous cases n+2 and 
n. 

n + 2 ones n ones 
A n +2 A n 
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This has sparsity 

sparsity of A n+2 + sparsity of A n + 2n + 2. 



dimension 2n + 1 = 2 k + 1 + 2i with i even: Here we combine the cases 
n = 2 k ~ 1 + 2% and 2 k ~ 1 + 1 and add a new top row. 



^2n+l — 



'2 fc-1 + 2t stars 2 fe ~ 1 + 1 stars' 



dimension 2n + 1 = 2 k + 1 + 2% with i odd: 



^2n+l — 



"2 fc_1 stars 2 fe_1 + 1 + 2i stars 



3. Unitary Matrices with First Row a Constant 

The researchers at the Frame Research Center needed unitary matrices 
with the first row a constant and maximal sparsity for their work on equian- 
gular fusion frames. They also wanted a large number of zero's in the matrix. 
For this algorithm, we can just alter the top row, then add a new top row 
and first column to the previous example. 

To make this more clear, we do two small examples and then the general 
case. First we do 3 x 3. 



l 



i -i 




V2 J 



Next we do 4 x 4. 




V2 V2 - 



Now for the general n x n case. 
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1 

n—1 



y/n(n-l) 



1 



1 



1 



sjn(n— 1) yjn(n— 1) ^Jn(n— 1) -y/ n(n— 1) ^/n(n— 1) 










/54 









/5-4 
3 






/5-4 
1 

/4^3 
2 



/5-4 
1 

/T3 
1 

1 



V5-4 
1 

1 

1 



4. Unitary Matrices with First Two Rows a Constant Modulus 

The researchers at the FRC needed these matrices for their work on con- 
structing tight fusion frames. 

Again, we give a few small cases to get this started. 

dimension 4: 



A A = 



-1 




1 


1 


1 


-1 


1 


1 


1 


1 


-1 


1 


1 




1 


-1 



dimension 6: 


















For Ag we can use a Hadamard again. But I was asked for "sparse" 
matrices doing this. 



dimension 8: 
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A* 











1 

8 

! 

12 










1 

8 

± 

12 









1 

8 

J_ 

12 

1 

6 

1 

2 

























1_ 

12 

2 
3 













12 

R 

6 
2 



And in general: 
A 2n = 



i_ 

2n_ 

_i_ 

2n 

71—1 



n(n— 1) 








n(n— 1) 



n-1 
n 



2n_ 

J_ 

2n 



1 

n(n— 1) 



5. Unitary Matrices: Made from 2 Constants 



2n 

J_ 

2n 



n(n— 1) 



We will look at various ways to construct unitary matrices from just two 
constants. 

5.1. First Case. In this section we construct a very special class of unitary 
of the form 



1 

n 



—b a a a ■ ■ ■ a 

a —b a a • • • a 

a a — b a • • • a 

a a a a • • • —b 
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For the case of R this looks like: 

1 
3 

For the case R 6 this looks like: 



-1 


2 2 


2 


-1 2 


2 


2 -1 



Or in a general form, the above looks like: 



-2 


1 


I— i 


1 


1 


1 


h- 1 


-2 


1 


1 


1 


1 


h- 1 


1 


-2 


1 


1 


1 


h- > 


1 


1 


-2 


1 


1 


h- > 


1 


1 


1 


-2 


1 


h- > 


1 


1 


1 


1 


-2 



2 
2 
2 
2 



2 
2 

-4 
2 
2 
2 



2 
2 
2 

-4 
2 
2 



2 
2 
2 
2 
2 

-4 



And for the general case of R n we have: 



1 

n 



-(n-2) 2 2 
2 -(n-2) 2 
2 2 -(n-2) 



2 
2 
2 



2 2 2 ••• -(n-2) 

5.2. Second Case. In this section we construct a class of unitary matrices 
of the form 



—a 


a 


a 


a 


-b 


b 


b 


b 


a 


—a 


a 


a 


b 


-b 


b 


b 


a 


a 


—a 


a 


b 


b 


-b 


b 


a 


a 


a 


—a 


b 


b 


b 


-b 


-b 


b 


b 


b 


a 


—a 


—a 


—a 


b 


-b 


b 


b 


—a 


a 


—a 


—a 


b 


b 


-b 


b 


—a 


—a 


a 


—a 


b 


b 


b 


-b 


—a 


—a 


—a 


a 



The rows and columns of this matrix are clearly orthogonal. To make the 
rows and columns square sum to one, we just need to have that 4a 2 +46 2 = 1. 

and a = 2b work. The above matrix 



20 



For example, in the above case, b 
represents a generalization of the Hadamard construction. In general it looks 
like: 
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Theorem 5.1. Let A,B be n x n unitary matrices. Choose a,b so that 
a 2 + b 2 = 1. Then the following matrix is unitary: 



C 



aA bB ' 
a A -bB 



Proof. A direct calculation gives that C*C = I. 

More general than this is: 
Theorem 5.2. Let A\,A%, . . . ,A n be m x m unitary matrices and let 

(. a ij)7j=i ' ^ e a un ^ ar y matrix. 
Then the following matrix is a nm x nm unitary matrix: 



□ 



B 



auAi ai 2 A 2 

CI21A1 022^2 



a-2nA n 



a n \A\ 02n^2 ■ ■ ■ 

Proof. A direct calculation gives that B*B = I 



□ 



Remark 5.3. The above theorem can be generalized to tight frames - which 
we do in a later section. 

5.3. Third Case. In this section we find unitary matrices of the form 



a 


a 


a 


-b 


-b 


b 


a 


a 


a 


-b 


b 


-b 


a 


a 


a 


b 


-b 


-b 


-b 


b 


b 


a 


a 


a 


b 


-b 


b 


a 


a 


a 


b 


b 


-b 


a 


a 


a 



For this to work, we need: 



3a 2 + 36 2 = 1, 



and for the inner products to be zero, 



3o 2 - b 2 = 0. 



Solving we get: 
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6. Two Tight Frames 

Here we will give a general construction for unit norm 2-tight frames. The 
FRC needed such examples made up of a single vector. Unit norm 2-tight 
frames are the row vectors of a 2n x n matrix having the following properties. 

(1) The rows square sum to one. 

(2) The columns square sum to two. 

(3) The column vectors are orthogonal. 

We start with a concrete example and then go to the general case. 



If we let 




We can think of reversing the rows of A to get the matrix: 

y/7 y/E y/3 -y/l 

V7 -y/E Vs Vl 
-y/7 \/5 V% VI _ 

We can also think of reversing the columns of A to get: 

' Vl V3 V3 -V7~ 

v 7 ! -\/3 VE V7 
-Vl v% VE V7 . 




A c = ^ 
16 
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If we reverse the rows of A and then reverse the columns of the resulting 
matrix we get: 




So our first matrix is of the form: 



Let us do this one more time. 



" A n 

A RC 





Vis 


VTi 


V9 


-V7 


VE 


Vs 


Vi 1 


VIE 




VTi 


V9 


V7 


-V5 


Vs 


VT 


VIE 


Vis 


-Vn 


V9 


V7 


VE 


-Vs 


VT 


VIE 




VTi 


-V9 


V7 


VE 


Vs 


-VT 


-VIE 




VTi 


V9 


V7 


~VE 


-Vs 


-VT 


VIE 


-Vis 


Vn 


V9 


-V7 


VE 


-Vs 


-VT 


VIE 


Vis 


-Vn 


V9 


-V7 


~VE 


Vs 


-VT 




Vis 


Vn 


-V9 


-V7 


-V5 


-Vs 


VT 


-VT 


Vs 


VE 


V7 


-V9 


Vn 




VIE 


Vi 


-Vs 


VE 


V7 


V9 


-Vn 




VIE 


Vi 


Vs 


-VE 


V7 


V9 


Vn 


-v^ 


VIE 


Vi 


Vs 


VE 


-V7 


V9 


Vn 


VlS 


-VTE 


-Vi 


Vs 


VE 


Vr 


V9 


-Vn 


-v^ 


-VTE 


Vi 


-Vs 


VE 


V7 


-V9 


Vn 


-Vis 


-VTE 


Vi 


Vs 


-V5 


V7 


-V9 


-Vn 


Vis 


-VTE 


. Vi 


Vs 


VE 


-V7 


-V9 


-Vn 


-Vis 


VIE _ 



In our earlier notation, this is of the form: 



A B 

A -B 

B RC j^RC 

B RC -A RC 



This construction will always work as long as the top row is an arithmetic 
progression. So the general case here looks like (in the 8x4 case): 
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— V a + 36 


Va + 26 


Va + 6 


Va 




y/a + 36 


-Va + 26 


Va + 6 


Va 




Va + 36 


Va + 26 


-Va + 6 


Va 


1 


Va + 36 


Va + 26 


Va + 6 


Va 


V4a + 66 


-\fa 


Va + 6 


Va + 26 


- Va + 36 




\fa 


-Va + 6 


Va + 26 


Va + 36 




yfa 


\J a + 6 


-Va + 26 


Va + 36 




yfa 


Va + 6 


Va + 26 


-Va + 36 



The columns are orthogonal, the rows square sum to 1 and the columns 
square sum to 2. i.e. This is a 2-tight frame. We can then iterate this 
procedure as above to get ever larger examples. 

7. Two Tight Frames: Weight in Front 

The idea now is to create 2-tight frames which have a big portion of 
their weight at the beginning. A group at the FRC wanted this for their 
work on the Kadison-Singer Problem. This will give an alternative concrete 
example of non-2-pavable projections. In [TJ it was first shown that the class 
of projections with constant diagonal 1/2 are not two-pavable. Later, in [2] 
a concrete construction of such projections was given using variations of the 
discrete Fourier transform. 

Let us start with an example and then see that it is a very general concept. 



"-Vl5 


VTi 


VT3 


VT2 


-Vi 


V3 


V2 


VT 1 


Vl5 


-VTi 


VT3 


VT2 


Vi 


-V3 


V2 


VT 


Vl5 


VTi 


-VT3 


VT2 


Vi 


V3 


-V2 


VT 


Vl5 


VTi 


VT3 


-VT2 


Vi 


V3 


V2 


-VT 


-Vl5 


VTi 


VT3 


VT2 


Vi 


-V3 


-V2 


-VT 


Vl5 


-VTi 


VT3 


VT2 


-Vi 


V3 


-V2 


-VT 


Vl5 


VTi 


-VT3 


VT2 


-Vi 


-V3 


V2 


-VT 


Vl5 


VTi 


VT3 


-VT2 


-Vi 


-V3 


-V2 


VT 


-VI 


V2 


V3 


V4 


-VT2 


VT3 


VTi 


VT5 


VT 


-V2 


V3 


Vi 


VT2 


-VT3 


VTi 


VT5 


VT 


V2 


-V3 


V4 


VT2 


VT3 


-VTi 


VT5 


VT 


V2 


V3 


-Vi 


VT2 


VT3 


VTi 


-VT5 


-VT 


V2 


V3 


Vi 


VT2 


-VT3 


-VTi 


-VT5 


VT 


-V2 


V3 


Vi 


-VT2 


VT3 


-VTi 


-VT5 


VT 


V2 


-V3 


Vi 


-VT2 


-VT3 


VTi 


-VT5 


. VT 


V2 


V3 


-Vi 


-VT2 


-VT3 


-VTi 


VT5 . 



We note that in the matrix above, the columns are orthogonal and the 
columns square sum to 2 and the rows square sum to 1. 

Now we will seriously generalize this example by understanding why it is 
working. What makes this work comes from the proof of how we get the 
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formula for 1 + 2 + 3 + • • • + n. To do this, we write these numbers in reverse 
and add them in pairs 

"1 2 3 ■ ■ ■ n - 2 n - 1 n 
n n-1 n-2 ••• 3 2 1 

We note that the sum of the columns are all the same and they add up 
to twice the sum of the rows. So we will work with integers n\ > ri2 > • • • > 
n 2 k-i and decide what we want the sums to equal, say m, and consider the 
matrix 



yjm — n\ yjm — n 2 



yjm — n 2 k-i 



yjm — n 2 k-i 

y/n 2 k-l 



sjm — n 2 k- 



^/n 2 k-x_ x 

Then we write the first row 2 k times and below it the second row 2 k 
times. We now need to assign signs of ± to this 2 k+1 x 2 k matrix so that the 
columns are orthogonal. This is easy to do. We use a variation of Hadamard 
matrices, i.e. 

"-11 1 1 

I- 111 

II- 11 

III- 1 



And by induction, 



A 



n+l 



An. At 



A r , 



~ Ar 



We use this twice, i.e. for each set of 2 k copies of the two rows. With 
this class of signs, we get a 2 k+1 x 2 k matrix with orthogonal columns and 
square row sums 

2&— l l 

(1) ^ra;+^(m-ra;) 



2 fc -V, 



i=l 



i=l 



and square sums of the columns is just 2 k m. So after normalization, this is 
a unit norm 2-tight frame. 



Our examples above do this. Also, the following works - for example: 



32 31 30 29 25 24 23 22 18 17 16 15 11 10 9 8 
8 9 10 11 15 16 17 18 22 23 24 25 29 30 31 32 
Further generalizations: 

Remark 7.1. Note that the terms may not all be decreasing above. This 
would require that 

n 2 k-i > m — n 2 k-i . 
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That is, 

m < 2n 2 k-\ . 

Remark 7.2. Note that we do not really need the terms to be decreasing 
above. All we really need is for the rows to sum to 2 k ~ 1 m and the columns 
to sum to m. For example, 

"20 24 19 25 1 7 2 6" 
6 2 7 1 25 14 24 20 

Remark 7.3. Note the we do not need the numbers to be distinct above. 
We can write one of them as many times as we like as long as we satisfy 
Equation [TJ 

Remark 7.4. Note that this gives another solution to the paper of Casazza, 
Fickus, Mixon and Tremain [2] where they used variations of the DFT to 
construct non-2-pavable matrices. A careful choice of the numbers above 
putting most of the weight in the first n — 1 positions then works in their 
proof. 

Here is another outline for constructing 2-tight frames. Fix a, b, c, al and 
let m = a + b + c + d. 

a b c d m — dm — cm — bm — a 

m — lm — bm — cm — d d c b a 

Then, the sum of the column numbers is a + b + c + d and we have 8 of 
these so the total sum after we turn this into a matrix is 8m while the sum 
of the rows is Am. So, after normalization, this is a unit norm 2-tight frame. 

8. General Tight Frames 
We will generalize one of our earlier examples. 
Theorem 8.1. Let {^4j}™ =1 be r x k{ matrices with orthogonal columns. Let 
( a ij)i=i jZi > ^ e a matrix with orthogonal columns. 

Then 

anAx ai 2 A 2 ■ ■ ■ a\ n A n 
ai\A\ a 22 A 2 ■ ■ ■ a 2n A n 

a miAi a m2 A2 ■ ■ ■ a mn A n 
is a mr x L matrix with orthogonal columns where 

n 

1=1 

Proof. Obvious. □ 
This leads to: 
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Theorem 8.2. Let A±, A 2 , . . . , A n be km x m matrices with unit norm rows 
and orthogonal columns which square sum to km. i.e. These are unit norm 
k -tight frames. Let 

( a ij)7j=i ' ^ e an orthonormal matrix. 

Then the following matrix is a kmn x mn unit norm tight frame. 

ai\A\ a\ 2 A 2 ■ ■ ■ a\ n A n 
a 2 \Ai a 22 A 2 ■ ■ ■ a 2n A n 

_a n \Ai a 2n A 2 ■ ■ ■ a nn A n _ 

Proof. Obvious. □ 

Remark 8.3. The same theorems hold with rows replaced by columns. 

9. Unbiased Bases 

Here is a way to get three unbiased bases in M 4 . One will be the unit 
vectors. The other two are below where +, — means ±1: 



And the last basis: 





+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 






+ 








+ 




+ 
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